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Abstract 


This paper presents an image method algorithm for the 
derivation of elastostatic solutions for point sources in bonded 
half spaces assuming the infinite space point source is known. 
Specific cases have been worked out and shown to coincide with well 
known solutions in the literature. 



Introduction 


Point sources (Green's functions are point sources) for some 
given (linear) governing differential equations and boundary 
conditions are important because of two main reasons. First, any 
localized process when viewed from a sufficient distance can be 
modelled as some suitably chosen point sources. Second, Green's 
functions can be used to reframe the governing differential 
equations and boundary conditions in an integral equation form; the 
integral equation form can, for example, be used as the basis for 
numerically analyzing a large class of problems using the boundary 
element method. 

This paper presents a new algorithm for the derivation of 
point sources of elastostatics in bonded halfspaces assuming the 
point source in infinite space is known. The method is similar to 
the image method that is familiar when deriving Green's functions 
in plane layered media where there is only one unknown scalar field 
in the governing equations such as in heat conduction, potential 
flow and electrostatics problems. In a sequel paper, the algorithm 
is then used to formally derive new Green's functions for any point 
source in a region consisting of an elastic layer perfectly bonded 
to two elastic half spaces. Numerical solutions for the displacement 
fields of nucleii of strain in an elastic plate are also presented 
in the sequel paper. 

Background 

There are many known Green's functions for half space problems 
in elastostatics. Most of the known Green's functions are 
specialized for a single half space having a stress free surface (a 
special case of bonded elastic halfspaces when one of the regions 
has zero rigidity) . Some of these known solutions are briefly 
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surveyed with occasional comments on the method of derivation. 

The most used point source solutions are the point force, the 
dislocation and the nuclei of strain (or double couple) solutions. 
The point force solution for 2-D plane problems in a half space with 
a free surface (Mellan 1932), 3-D problem in a halfspace with a 
free surface (Mindlin 1936) and 3-D problem in bonded elastic 
half spaces (Rongved 1955) are known. Rongved obtained the Green's 
function through the use of the Papkovi ch-Neuber potentials and 
arguments from harmonic analysis; the resulting solution is in the 
form of the sum of a point force solution in infinite space and 
some point sources at the image point with respect to the interface 
plane. 

The elastic fields of screw and edge dislocations in bonded 
halfspaces were first given by Head (1953 a,b) . The screw 
dislocation problem is obtained by the method of images (since 
there is only one field variable) . 

There are six nucleii of strain sources. The solution to the 
first (double couple in a plane parallel to the free surface) was 
given by Steketee (1958), the remaining five sources were given by 
Maruyama (1964). Maruyama used image nucleii of strain sources to 
cancel the tangential component of the surface traction on the free 
surface. He then used the Boussinesq solution (in Galerkin vector 
representation) and the remaining normal tractions on the free 
surface in a Hankel/Fourier transformed space to obtain the rest of 
the fields afterwhich he transformed the solution back to real 
space. This procedure is highly specific to half space problems 
with a free surface and cannot be generalized to multiple layered 
systems. 


Finally, we note the existence of an image method for 
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perfectly bonded elastic halfspaces in terms of the 
Papkovi ch-Neuber potentials (Aderogba 1977). Aderogba presented the 
algorithm for obtaining the four image potentials which involves 
multiple integrations with respect to the coordinate perpendicular 
to the interface plane and differentiation with respect to all 
three coordinates. The algorithm based on the Hansen potentials 
presented in this paper involves 3 potentials only, and only 
differentiation of the (infinite space) potentials with respect to 
the coordinate perpendicular to the interface plane (as well as 
multiplication by scalars) is required to obtain the image 
potentials. This distinction is especially important when the image 
algorithm is repeatedly applied to obtain the fields due to point 
sources in regions consisting of an elastic layer perfectly bonded 
to two elastic halfspaces. 

Preliminary Considerations 

The image method presented in this paper is dependent upon 
expressing the displacements in terms of potentials. The specific 
potentials employed are the analogue to Hansen's potentials for 
elastostatics and dynamics. Unlike Ben-Menahem and Singh (1968) the 
potentials are not expanded in terms of eigenfunctions; instead the 
algorithm operates directly on the potentials. Note however, that 
the eigenfunction expansion technique was used in the derivation of 
the algorithm (see Appendix 2). 

Specifically, we express the displacement field is expressed 
in terms of the Hansen potentials Y and ‘ t * ie following 

manner: 


u(h,6,£ R ,£ L ) = Nfh,^) + F(6,h,y> 2 ) + M(h,Y> 3 ) 
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N(h,y , 1 ) = v?> 1 (x,y,z-h) 

E(o,h,r 2 ) = + 2-e 2 -|^ 2 (x,y,z-h) 

- vr 2 (x,y,z-h) 

- 2*5* (z-h) r|j(’ 2 (x,y,z-h) 


( 1 ) 


£1(h / ¥ > 3 ) = v x [e^ 3 (x,y,z-h) ] 


where: 


v is the gradient operator 

v x is the curl operator 


v 2 y > 1 = v z ? 2 = v i r > 3 = 0 
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A 


_ A + fj 
~ T~+~3jj 

is the Lame constant 


fj is the shear modulus 

h is a scalar for shifting the z-coordinate 


Note that the potentials * 1# ? 2 and ? 3 have to be harmonic in 

order for N, F and M to satisfy equilibrium. The Cartesian 

components for the displacements, strains and stresses are given in 

appendix 1 of Fares and Li (1986). It is shown in appendix 3 that 

p is associated with the antiplane mode of deformation. 

L 

In order for these potentials to be useful, a method to obtain 
these potentials given an elastic field satisfying equilibrium is 
described below. Note that: 
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> 

V-N = v-M = 0 V-F = 2- ( 1-5 ) * 

az 2 

v x N = 0 

( 2 ) 

d 2 V . d 2 ! 9 - 

vxE = 2- (!+«)• -*e - 2- (1+5) 

dydz X axdz ^ 

a^ 3 . a 2 ^ . a^ 3 . 

v x M = • e + • e + • e 

axaz x ayaz ^ az 2 z 


Hence for a given displacement field u, the following may be 
calculated: 


V-U = 2 - (1-6) ' 


a 2 y>. 


az J 


and 


(3) 


d 2 r. 


(v x u) -e = 


az 2 


(4) 


The potentials tp and r may be obtained by integrating (3) 


and (4) : 


*2 = 


vu 


dz dz + z-F (x,y) + G (x,y) 

J 2 - ( 1 - 6 ) 


(5) 


y 3 = 


dz 


dz[ (v x u) *e z + z-F 3 (x,y) + G 3 (x,y) 


The integration constants F^'s and G^'s are chosen such that 
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^2 and p are harmonic in the required region. In addition, all 
singularities of the potentials must be in the region where the 
source occurs. This is made clearer in appendix 3 when we consider 
examples of the use of the algorithm. Finally, once and are 
determined, whatever remains in the displacement field (see 
equation 1) is ascribed to If the given displacement field does 
satisfy equilibrium, the field should be expressible in terms of 
these three potentials (see Ben-Menahem and Singh 1968, and Morse 
and Feshbach 1953 ) . 

The Hansen potentials for a point force, and a line force 
perpendicular to the z-direction are given in Appendix 1. 

The Algorithm 

The algorithm and the notation associated with it can now be 
described. Consider two elastic half spaces perfectly bonded along 
an interface plane at z=0 (see figure 1). The material properties 
of region 1 are described by fi and 5^ , and of region 2 by and 
• Next we define the following: 


ip(x,y,z) = TP(x,y,-z) 


s Ns/Ph 


( 6 ) 


= llilll 

TK+tT 


b = 


(6,+l) 

(t6 2 +1) 


Note that if f is harmonic then T is also harmonic and hence 
can be used as a Hansen potential for N, F and M. 

The algorithm states that if we have the representation for a 
point source in infinite space of elastic constants similar to 
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those of region 1 at the location x=y=0 and z=h described by the 
displacement field: 


u° = B°(h,a lf £°,E°) 


( 7 ) 


then the displacement fields in regions 1 and 2 for a similar point 
source in region 1 at x=y=0 and z=h are given by: 


u 1 = u° + uf-h,^, £*,£*) 
u 2 = u(h,6 2 ,£ 2 ,n£) 


( 8 ) 


where the image potentials may be obtained from the source 
potentials with the following operations: 


E* = R R (-h,a,b,s 1 )-y° 






^ R (-h,a.,b,5 1 ) 



-Zs.d-alh-lj 

n O 

+ ( l-b) - 46 f ( l-a)h 5- 

dz 


+ (l-a) 

+ 26 x (l-a)h-| i - 


+a 


- 2 ( 6 2 



b-6 i a ) h ’lz 


+b 


| R (h,a,b,6 2 ,6 x ) = 
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1l {7 > 


1-7 

ITT 


|.( 7 ) 


1 

+2 

1+7 


( 9 ) 


The above operators are denoted by 'R' and 'T' and stand for 
'Reflection' and 'Transmission' operators respectively, in analogy 
with wave reflection and transmission operators for plane waves in 
elastodynamic problems. 


We note that: 

= Z R (-h,a,b,6 1 )-Y^ = ^(-h,a,b,6 l ) 
and: (10) 

^(■•h.a.b.c,) = R R (+h,a,b,6 l ) 

12 0 

Note that the £_ ' are simple multiplicatives of £ . The case 
L L 

when £^ = 0 corresponds to the purely anti-plane problem, and thus, 

the algorithm reduces to the scalar image method for that case. 

The derivation of the above algorithm is given in Appendix 2, 
and some sample known solutions are rederived in appendix 3; namely 
the screw dislocation in a half space with a free surface and 
Mindlin's solution of a point force interior to a half space. 

Conclusions and further recommendations 


A vector image method has been presented, for elastic problems 
with planar interfaces. An algorithm has been presented on how to 
derive point source solutions for two bonded elastic half spaces. 
Specific cases have been worked out (in appendix 3) and shown to 
coincide with well known solutions in the literature. Further 
details and sample cases could be found in Fares and Li (1986) . 
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The method of deriving the algorithm suggests that an 
analogous algorithm can be obtained for spherical interface 
problems, and 2-D (but not 3-D) cylindrical interface problems in 
elastostatics . This suggestion is supported by the existence of a 
scalar image method and Hansen potential representations for both 
these geometries. 

Finally, it would also be of interest to investigate 
equivalent algorithms for other governing equations. For example, 
elastodynamics and poroelasticity could be potential candidates for 
such an investigation. Elastodynamic problems, in particular, do 
have Hansen potential representations that have been well 
established and used and could be investigated first without the 
considerable preliminary formulations that are needed for 
poroelastic problems. 
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Appendix 1 : Sample potentials for some point sources 
Point Force: 


The displacement field due to a point force at the origin can 
be written as: 


u. = 


1 4np~’ (1+5) ' 


P i4- + 5 -Pk- X k' X i- ^3 

r 


(i.i) 


where: 


6 s A+/J 


A +3 /J 


1 

K 


3-4u 


2 2 2 2 
r = x + y + z 


p,. are the magnitudes of the point forces in the 

. th , . . 

a -direction 


The Hansen potentials for the point force can be obtained by 
using (1) and (5) : 




Pirlz + p 2'r5z * P 3 ' ln(r±2) 




Pa'iS 1 p 3' ln(r±z > 


( 1 . 2 ) 


1P 3 =p- 


p l* ( 1+5 )'riz ” p 2* ( 1+5 ) 'rlz 


where: 


P = l/[4»/u* ( 1+6 ) ] 


Note that if the upper (lower) "sign" of ± is chosen in one 
expression, the upper (lower) "signs" must be chosen throughout for 
all the potentials. Also note that taking r+z (r-z) in the 
expressions makes the potentials (but not necessarily the 


displacements) singular when x=y=0 and z<0 (z>0) . 


Line forces at x = 0 acting on the z = 0 plane 


The displacement field due to a line force can be written (for 
plain strain) as: 


u i 4np5 


-p.-lnc + «'V 3! k- x i- 


and 


for i,k = 1,3 


(1.3) 


u 


2 


= 0 


where : 


A+u 

a = — 

A +2 fj 


5 s ^ 0 - 
A+3p 


C 


2 



+ 


„2 


p and p are the magnitude of the line forces 

1 3 

The Hansen potentials for the line force can be obtained by 
using ( 1 ) and ( 5 ) : 


r 


l 


a 


a 

SnJjE 


~ Pi 


"Pi 


- Pi 


z-arctan(— ) - x-ln? + (1+5) -x 
z- In? - z + x-arctan(§) 

z-arctan(9) - x-lnC + (1+5) -x 

X 

z- In? - z + x-arctan(^) 
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Dislocations parallel to the 2=0 plane: 

The displacement due to a dislocation along the y-axis (plane 
strain) can be written as: 

r ^ ■ 

u i * [ sr arctan( I> ' ‘ik'sr ln{ _ 

+ ‘ <2f,t ik d k ) ' :ln5 + . 

for i,k,n = 1,3 

and: (1.5) 


where : 


'ik 


+1 for i = 1 , k = 3 

-1 for i = 3, k = 1 

0 otherwise 


d 1 and d 3 


are the slip magnitude of the dislocations 


We note that the terms in the second brackets expressing the 
displacements are of the form of line force expressions with 
equivalent magnitudes of 2//^^ and thus their Hansen potentials 
are already known. The Hansen potentials for the terms in the first 
bracket can be shown to be: 



z-lnt - 2 + x-arctan(^) 


+ d 


3' 


2 

2 -arctan(— ) 
x 


x-lnf + x 
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^ioracicet _ 


JMet 


( 1 . 6 ) 
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Appendix 2: Derivation of the vector image method algorithm 

The eigenfunction expansion method for elasticity problems in 
layered media was first formulated by Ben-Menahem and Singh 1968. 

We have vised this method to derive the algorithm discussed in this 
paper. The notation (as far as possible) is the same as in the 1968 
reference paper, although some new temporary terms have been 
defined in order to simplify the algebra for this specific 
implementation . 

Any elastic displacement field satisfying the equilibrium 
equations : 


v J u + (1+A/jU) -vvu = 0 (2.1) 

can be written as the sum of N, F and M (see (1) ) . 

Using the method of the separation of variables in cylindrical 
coordinates on the potentials y 5 , , and ^3 the form: 

f = Z(z)-R(r)-F(e) 

We get: (2.2) 

V = exp(±kz) • J m (kr) •ejp(±im 6 ) 

where J is the Bessel's function of the first kind and of m 
m 

order. 

Reexpressing r , r and in the above form and carrying out 

X £ 

the v and r x and X- operations in ( 1 ) , we get : 

oZ 



m 



A -N 
m m 


+ B +, F + 
m m 


+ B'-F + C + -M + + C -M 
mm mm mm 


•dk 

(2.3) 
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where: 


A*, B ± and C* are constant coefficients of U, F and 
m m m 

M dependent on 'm' only. 


N* = exp ( ±kz ) • 

Tl 


±P + B 
“m “m 


= exp(±kz) • (±l-28kz) - (l±26kz) -B 


-m 


and: 


= exp ( ±kz ) • C 


~m 


P = e • J ( kr ) • exp ( ime ) 
~m z m 


B = (e -J?— + ^-)J (kr) -exp(ime ; 

-m r c»kr e kr ae ' m 


C = (e • - e )J (kr) -exp( ime ) 

-m ' r kr ae e akr m 

(2.4) 

In the above expressions for P. B and C there is the 

-m mi Tti 

implicit understanding that we can consider either the real or 
imaginary components of the expressions seperately. 


From the above expressions for the displacements, the 
expressions for the tractions at a plane z=constant can be found, 
and we rewrite the above as: 


oo 

- = Z 

m=0- 

00 

f-l 

m=0 


(u? + uj!) -dk 


(T* + A-dk 
mn mi 


(2.5) 


(2.6) 


where: 
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and: 


u = x *P + y -B 
-m m -m 2 m “m 




= 2k-X_*E m + 2k-Y m -g m 


m 


u = z -C 
-m m -m 

i = k-Z*C 
ti m 


(2.7) 


x m = A**exp(kz) - A ffi *exp(-kz) 

+ B+- (l-25kz) -exp(kz) + B • (-l-25kz) • exp(-kz) 

y m = A*-exp(kz) + A~-exp(-kz) 

+ B** (-l-25kz) -exp(kz) + B~- (-l+25kz) -exp(-kz) 

z = C + -exp(kz) + C -exp(-kz) 
m m ^ m c 


X = A'-u-exp(kz) + A -/u-exp(-kz) 
m m m 

+ B + -u6- (l-2kz) -exp(kz) + (l+2kz) -exp(-kz) 

m m 


Y m = A m -/J-exp(kz) - A m -/j-exp(-kz) 

+ B + -u6- (-l-2kz) -exp(kz) +B -p5 • ( l-2kz) • exp(-kz) 
m m 


Z m = c m -A /-exp ^ kz ) “ C m -/u-exp(-kz) 

( 2 . 8 ) 

R L 

Notice that the components are uncoupled from the In the 

sense that the A ± / B* coefficients do not affect the components 
m m ~ 

and the C coefficients do not affect the u" components. We 
m ~ 

R L 

therefore treat the u and the u components seperately when 

■“iTl “HI 

analyzing a specific problem in terms of the Hansen potentials. 
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Consider the specific geometry shown in figure 2.1. The region 
consists of two elastic materials seperated by a planar interface. 
The material elasticity parameters used to characterize the regions 
are taken to be /j x , S x and 6 2 . A point source exists at the 
position z=-h. The problem is to find the displacement fields for 
region 1 (z<0) and for region 2 (z>0) under the influence of the 
point source, such that the displacements and the tractions are 
continuous across the interface plane (z=0). 

R L 

In what follows, we are manipulating u and u in equation 

“m “th 

2.5 for a fixed 'm' , but the 'm' subscript will be dropped for 
brevity. First, we express the displacement and traction (on a 
z -plane) for (z+h) > 0 (which includes z=0) of a point source of 
arbitrary nature (using the eigenfunction expansion method and 
expressing the m'th component in matrix form) in the following way: 


U(P) 


r -i 

1 

-1 - 26 x k-(z+h) 


1 

0 < 

u(B) 


+i 

-1 + 25,k- (z+h) 


B°- 

T(P) 


+2k,u x 

+2ty x 6 x • [l+2k- (z+h)] 



T(B) 


-2k^i 

+2k,u x 5 x • [l-2k- (z+h) ] 



JO L 


u(C) 


+1 

C° *e? 5 )(-k|z+h| ) 

T(C) 

0 




(2.9) 


and we define: 
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_fl_ 


u(P) 


JfL 

= —1 . 

u(B) 



— J. * 

T(P) 


. S 4 . 


T(B) 

0 


[•.1 

— * . 

u(C) 


_ S 5 . 

— X * 

T(C) 




L 

'o 


( 2 . 10 ) 


Now the elastic fields in region 1 are expressible as: 


U(P) 


+1 

. 

+1 - 26 2 


A 1+ ' 

u(B) 


+1 

-1 - 26 1 k , z 


B 1+ 

T(P) 1 



+2k/i 1 6 1 • (l-2kz) 

1 

. 

T(B) 

1 

+2k^ ; 

. 

+2k^/ 1 6 1 • (-l-2kz) 



•exp(kz) 


+ 


terms due to the ] 
point source j 
as given above j 


u(C) 


+1 

T(C) 

1 

+1 gj 


•exp(kz) 


+ 


terms due to the 
point source 
as given above 


( 2 . 11 ) 


And the elastic fields in region 2 are expressible as: 


u(P) 


-1 

-1 - 26 2 kz 

u(B) 


+1 

-1 + 26 2 kz 

T(P) 


+2ty z 

+2k,u 2 5 2 • (l+2kz) 

T(B) 

i 

N 

£ 

CM 

1 

+2kfi 2 5 z • (l-2kz) 


J 2 L 




•exp(-kz) 
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ORIGINAL PAGE IS 

Of POOR QUALITY 


U(C) 


+1 


'c 2 -' 

T(C) 

2 

-k^ 2 




exp(-kz) 


( 2 . 12 ) 


Applying the condition that and are to be continuous 
(for each m) along the interface plane z=0, we get: 


+1 

+1 

+1 

+1 


’ A 1+ ‘ 


S 1 ’ 

+1 

JL 

-1 

+1 


B 1+ 


S 2 

+2k^ 

+2k ) u 1 6 1 

N 

CM 

1 

-2k^ 2 6 2 


> 1 
roi 
l 


S 3 

+2k^x 

-2^i «i 

+2kv 2 

-2k|U 2 6 2 


B 2 

• 


U) 1 
4* 


f« 

■ 


' c 1+ ‘ 


1 

if) 

(J1 

1 

i 

+ 

£ 

j 

* 

c 2 " 

j 

1 

6 


1 + 1 + 1 + 2 - 2 - 2 - 

Now we solve for A , B , C and A , B , C by inverting 
the 4x4 and 2x2 system of equations. We obtain: 


i 

> 

+ 


B 1+ 

1 

1 

CM 

c 

A 

r ~ — — 

Cd 

ro 

l 


’ c 1+ ‘ 

i 

o 

ro 

1 

7+1 ’ 

L 



+d/ 2-7 -5 j 

+^/ 2 - 7 - 6 j 

7+6, 

21 Vjl 

7+6, 

2k^ 


V 

+-d/ 2— 75 A 5 2 —6 A 

— ^/ 2+76 x 6 2 +6 j 

76 2 +l 

2^1" 

_76 2 +l 

2i^~ 


CM 

U) 

+7 5 L 5 2 +5 L 

-75 1 6 2 -6 l 

76 2 +l 

_2iz. 

76 2 +l 

!ieez_ 


s 3 

+^+Si 

+7+6, 

1 

7+5 : 
2tyi 

7+6, 

2k^, J 


. S 4 . 


+7 

+ K7 


V 



-1 


! 

. S 6 . 


(2.14) 
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where: i = 

A = 2*(7+6,)-{76 2 +l) 

Expressing the S's in terms of A° , and , and 
simplifying the expressions we get: 


V 


0 

3 

- b 



B 1+ 


1 - a 

26 i • ( 1- 

-a) -kh 


A 2 ' 


a 

25,a-kh 


B 2 - 

- J 


0 


b 



1 

o 

I-* 

+ 

! 

1 - 7 ' 

1 + 7 


'c 0 - 1 

•exp(-kh) 

c 2 - 

J 

2 

1 + 7 


L j 




•exp(-kh) 


(2.15) 


where : 


a = (S^D/b+Sj 


b a (6 l +l)/(7-6,+l) 


(2.16) 


Therefore we find that the displacement field (for a given m) 
in region 1 and region 2 can be written as: 


u(P) 


+1 

+l-26 1 kz 

u(B) 


r 1 

-1-26 ,kz 

j 

1 

L 




i-a 


1-b 
25 , • 


( 1-a) -kh 



•exp[k(z-h) ] 


+ 


source terms 
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u(C) 

= +1 


1 

U(P) 

-1 

u(B) 

+1 

J 

2 

u(C) 

= +1 


n 

L. *- 


*exp[k(z-h)J + source terms ] 


(2.17) 



+a 


0 


•exp[-k(z+h) ] 


;) = +1 • • c •exp[-k(z+h) ] 

(2.18) 

The u(C) terms are in a form from which we can deduce the 


i 

the 

now 

try 


+1 


+1 

a 

[+1 

az 

+1 

,2 

+1 

9 

dZ~ 

+1 


+1 +1-26 x k- (z-h) 


— -- -exp[k(z-h)] 


• exp [k( z-h) ] 


A c -exp[k(z-h) ] 


source terms 


u(P) 


-i 

-l-26 a k- (z+h) 

u(B) 

2 

+i 

-l+2S 2 k- (z+h) 


— -exp[-k(z+h) ] 
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-1 

-l-26 2 k- (z+h) 

+1 

-l+26 2 k- (z+h) 


®b 


•exp[-k(z+h) ] 


+ 


a 

dZ 


-1 
+ 1 


•exp[~k(z+h) ] 


( 2 . 20 ) 


Noting that: 
a n 

— — exp[k(z-h) ] = k n - exp[k(z-h) ] 
az 

and 

.n 

— ^exp[-k(z+h) ] = (-k) n -exp[-k(z+h) ] 
az 


( 2 . 21 ) 


We obtain: 


A + = (1-b) -B° B + = ( 1-a) -A° 

cl 3. 

aJ = -26, • (1-a) -h-A 0- + 46, • (1-a) -h-B 0- 

B* = 25, • (1-a) -h-B 0- A* = -46 , • ( 1-a) -h 2 -B°' 

and: 

A" = a- A 0- B“ = b-B°" 

a a 

A^ = 2 • ( 5 2 b-6 , a ) • h • B° _ B^ = A~ = 0 


( 2 . 22 ) 


(2.23) 


Since the above relations are true for each component of a 
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potential, then they must be true for the whole potential and we 
get: 


if: 


U° =N(-h,P°) +F(6 1 ,-h,?°) +M(-h.y°) 


then: 


a = 


u° + N(h, (1-b) -v^) + |^(h,-25, • (l-a) -h-v^) 

+ | 2 K(h,+4fi?.(l-a)*h-^) 

2 

+ — *N(h,-46 l • (l-a) •h^*F?) 
dz 2 

+ F(6 lf h,(l-a)^) 

+ | i F(5 l ,h,25 l -(l-a)-h-r^) 


+ M(h,^.y 0 


— -y°) 
1+7 3 ] 


u = 


N(-h,a^°) + ^(-h,2-(e a b-6 1 a)-h-?’°) 

«L o a-, £ 

+ F(5 2 ,-h,b?°) 

2 0 
+ M(-h >T i--^) 

'1+7 3' 


Noting that: 


(2.24) 


(2.25) 


N(h,cst-¥ > ) = N(h,cst- ?) 


az 


az 


n 


and 


l^’(6,h,cst-Y > ) = N(h, -25j • cst-^-5 0 ) + F(6,h,cst--^- 

d*- dZ dZ 


(2.26) 


where : 


cst 


is a constant 
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We obtain the algorithm given in the main body of this paper 
(two minor differences are: i) The statement of the algorithm in 
the paper considers region 1 to be at z>0 and hence z=+h instead of 
z=-h to be the location of the source point and ii) A formalism in 
terms of matrix operators is implemented in the main text ) . 
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Appendix 3: Derivation of some sample Green's functions 
through the use of the image algorithm 

In this section, we consider displacement fields in Cartesian 
components for some sample point source problems in a half space 
with a free surface. The solutions that will be rederived are 
readily available (and established) in the literature and hence 
serve as an empirical check of the algorithm. In addition, these 
specific examples help clarify details of the application of the 
algorithm. From the algorithm presented in the main text and from 
the examples to follow, it should be clear that many established 
elastostatic solutions and new solutions (in part II of this paper) 
can be obtained in a systematic fashion. 


By considering a halfspace problem with a free surface, we 
obtain the following simplifications: 


y = 0 


1-a = -l/fii 


call: 


6 = 6 , 


then: 


| R (h / a,b,6 l ) = 


i-b = - 6 ! 


• 

„ 2 1 


a 

-2h*- — 
oz 

-6 + 45h I. 2 - 





Rr(y) = 

1 

-1/6 

+2h 'lz 



. 

J 



(3.1) 

The following example problems will be considered: 

I. Screw dislocation '(Antiplane problem). 

II. Point force acting interior to the halfspace with a free 
surface (Mindlin's solution). 

i) The point force is in the x direction. 

ii) The point force is in the z direction. 
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I. Screw dislocation (antiplane problem) 


For the antiplane problem, all that is required is to obtain 
the image potential with respect to the interface plane since the 
matrix is the identity operator. Also, we notice that getting 
the image of a given M type (see main text, particularly equation 
8, and with specialization in 3.1) displacement field is equal to 
the M displacement field of the image potential describing that 

field (i.e. M(h,£) = £l(-h,£)), and hence we can directly operate on 
a given displacement field when using the algorithm for a purely 
antiplane problem. This corresponds to the scalar field image 
method for the antiplane case. 

As an example we consider the field due to a screw dislocation 
in the plane perpendicular to the x-z plane at location z=h and 
~0. The field due to the dislocation in infinite space is: 

u = arctan[ ( z-h ) /x] (3.2) 

The image field will be u which implies that the combined 
fields give: 


u = arctan[ (z-h)/x] - arctan[ (z+h) /x] (3.3) 

Of course this is but a simple application of the scalar image 
method. 

II . Point force acting interior to a half space with a free surface. 

For this problem the potentials for the point source in 
infinite space are given in appendix 1. However, we have a choice 
of where to locate the singularities of the potentials. Since we do 
not want the image potentials to introduce any new sources inside 
the half space, we choose the infinite space potentials to have all 
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their singularities in that halfspace. The following functions will 
need to be calculated: 


'a a a 


¥° 

[a a a a 2 a 2 a 2 ' 



ax' ay' az 


2 

. - 

ax' ay' az' axaz' ayaz' 2 

L o Z J 


i 


d 

a a a 


r_oi 

fi\ 


a a a 



az * | 

ax' ay' azj 


az 2 

ax' ay' az 

M 

L 2j 


a 

a a a a 2 a 2 a 2 ' 



az 

ax' ay' az' axaz' ayaz' 2 

L dZ J 


2 . 


There are some repetition in the suggested functions to be 
calculated since the partial differentiation operations are 
commutative when the function is sufficiently smooth. 

We note that the potentials to be differentiated are linear 
combinations of the following functions: 


? A = x/[2- (r+z) ] 
¥ q = -[ln(r+z)]/2 


( 3 . 4 ) 


In addition, we have an antiplane potential for this case 
which is a linear multiple of the following function: 


= y/[2- (r+z) ] 

and we will also have to calculate: 


d 1 


\p 

_ax' ayj 


B J 


( 3 . 5 ) 


Now we can perform the required differentiations and we also 
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note the following: 

(x 2 z)/[r 3 * (r+z) ] + (x 2 z)/[r 2 * (r+z) 2 ] = x 2 /r 3 - x 2 /[r* (r+z) 2 ] 
z/[r- (r+z)] = 1/r - l/(r+z) 

y 2 /[r- (r+z) 2 ] = 2/(r+z) -1/r - x 2 /[r- (r+z) 2 ] 

(xyz)/[r 3 - (r+z)] + (xyz)/[r 2 - (r+z) 2 ) ] = (xy)/r 3 - (xy)/[r- (r+z) 2 ] 
“5/2+1/ (25) = fJ’ (l+6)/(A+^) 

* (3.6) 

i) Case when the force is acting in the x direction: 

Defining: 

2 2 2 . . .2 
r 2 = x + y + (z+h) 

We get after simplification and the use of identities similar 
to those given in (3.6), but with z replaced by (z+h) and "r" 
replaced by "r 0 " wherever they occur: 

0 

u = u 

+ [l/4n/j(l+6) ] • 

e x - [5/2 + 1/(25) + 1 - 2- (1+5) + ( 1+5 ) ]/(r 2 +z+h) 

+ [-5/2 - 1/(25) - 1 + (1+6)] -x 2 /[r 2 - (r 2 +z+h) 2 ] 

+ [-1 + (l+6)]/r 2 

+ x 2 /r 3 + 26hz- [ 1/r 3 - 3-x 2 /r 2 ] i 
+ e • [-5/2 -1/(25) - 1 + (1+5)]- (xy)/[r 2 -(r 2 +z+h) 2 ] 

+ (xy)/r 3 + 26hz- [-3-x-y/r 3 ] | 

+ e z * [-5/2 +1/(26) ] -x/[r 2 * (r 2 +z+h) ] 

+ (z-h) -x/r 2 - 65hz-x- (z+h) /r 3 


(3.7) 
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Noting that: 

G/[4»ff* ( 1+6) ] = l/[16nfJ‘ { 1-v) ] 

1/6 = 3-4P (3.8) 

-6/2+1/ ( 26 ) = 46 • ( 1-v ) • ( l-2fc» ) 

we find that the above result (3.7) coincides with the solution 
first obtained by Mindlin (1936) and shown in Mura (1982). 

ii) Case when the force is acting in the z direction: 

After simplifications we get: 



+ [ 1/4«ai (1+6)]- 


e • 
x 


+ e 


+ e • 

z 


+ [5/2 - 1/(26) ] -x/[r 2 - (r 2 +z+h) ] 

3 5 

+ (z-h)-x/r 2 + 66hz-x- (z+h)/r 2 

[6/2 -1/(26) ]-y/[r 2 -(r 2 +z+h)] 

+ (z-h)-y/r 2 + 66hz-y- (z+h)/r 2 ] 

[6/2 +l/(26)]/r 2 

+ [(z+h) 2 - 26hz]/r 2 + 65hz- (z+h) 2 /r 2 


3.9) 


Noting the relations given in (3.8) and: 


6/2+1/ ( 25 ) = 6- [8- (1-u) 2 - ( 3-4v ) ] 


(3.10) 


we find that the above result (3.9) coincides with the solution 
first obtained by Mindlin (1936) and shown in Mura (1982). 
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Figure Cast ions 

figure 1: 

2 bonded elastic 

figure 2.1: 

2 bonded elastic 


halfspaces with a point source at z=h. 
half spaces with a point source at z=-h. 
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